Mechanism of magnetism in stacked nanographite 
with open shell electrons 
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Abstract 

Antiferromagnetism in stacked nanographite is investigated with using the Hubbard-type mod- 



> 

. els. The A-B stacking or the stacking near to that of A-B type is favorable for the hexagonal 

00 
O 



nanographite with zigzag edges, in order that magnetism appears. Next, we find that the open 
shell electronic structure can be an origin of the decreasing magnetic moment with the decrease 
of the inter-graphene distance, as experiments on adsorption of molecules suggest. 
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1. Introduction 



Nanographite systems are composed of the structural units with stacking of graphene sheets of 
the nanometer size [1]. They show novel magnetic properties, such as, spin-glass like behaviors 
[2] , and the change of ESR hue widths while gas adsorptions [3] . Recently, it has been found 
[4,5] that magnetic moments decrease with the decrease of the interlayer distance while water 
molecules are attached physically. 

In the previous papers [6,7], we have considered the stacking effects of hexagonal nanographite 
layers in order to investigate mechanisms of antiferromagnetism using the Hubbard-type model 
with the interlayer hopping integrals and the onsite repulsion U. We have taken account of the 
nanographite systems where the numbers of the sites and electrons in a layer are same, and 
therefore the electronic systems of the layer have closed shell structures. In the calculations, 
the finite magnetization develops, as the hopping interactions between layers increase in the 
case of the A-B stacking [6] . The same conclusion has been obtained when the system is near 
the A-B stacking and the interlayer distance becomes shorter [7]. The A-B stacking should 
exist in nanographite systems, because the exotic magnetisms have been observed in recent 
experiments [1-3]. The decrease of the interlayer distance while attachment of water molecules 
makes ti larger. However, it is known that the magnetism decreases while the attachment of 
molecules [4,5]. The calculation for the closed electron systems cannot explain the experiments 
even qualitatively. 

The purpose of this paper is to extend the previous calculations to the systems with open 
shell electronic structures in the isolated nanographite layer. There are several candidates for 
the open shell of electronic structures. One of them is the electron number changes due to the 
presence of the active side groups. The electron or hole is donated from the side group. The 
effects can be modeled by a site potential in the idea of the model hamiltonians. This idea will 
be investigated in the present paper. The other origin of the open shell is a geometrical origin. 
The phenalenyl molecule C13 with three hexagonal rings [8,9], and the triangulene C22 with 
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six hexagons [10,11] are the examples of small graphene layers which have open shell electronic 
structures. The stacking effects of such kinds of molecules will be reported elsewhere [12]. 

We shall study with the two kinds of models which have been used in the papers [6,7]. 
Site potentials which simulate the additional side groups are introduced in the Hubbard-type 
hamiltonians. The change of the electron number is taken into account, too. We will show that 
the A-B stacking or the relative interlayer relation near the A-B stacking is favorable for the 
hexagonal nanographite with zigzag edges, in order that magnetism appears. We also find that 
the open shell electronic structure can be an origin of the decreasing magnetic moment with 
adsorption of molecules. 

In the next section, we explain our models and review the method of the numerical calcula- 
tions. In sections 3 and 4, the results of the two models are reported and discussed. The paper 
is closed with a summary in section 5. 

2. Models 

First, we will investigate with the model I [6], 

H ^ -t y: E(cU>+h.c.) 

(i,j):intralayer 

(i j):interlaycr 

+ ^IZ^i.T^U' (1) 

i 

where rii^a = cl^crCi,a for cr =1 and j; Cj^o- is an annihilation operator of an electron at the ith 
site with spin a; the sum of the first line is taken over the nearest neighbor pairs {i,j) in 
a single layer of the nanographite; the sum of the second line is taken over sites where the 
distance between two positions of the neighboring layers is shortest; ti is the strength of the 
weak hopping interaction between neighboring layers; the positions of are shown by the filled 
circles in Fig. 1 (a); and the last term of the hamiltonian is the strong onsite repulsion with 
the strength U. We will vary the strength in a reasonable range. The actual values of U for 
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carbon atoms are of the similar order of magnitudes as in this paper: for example, we have 
found U = At for the neutral Cgo and C70 [13] and U = 2t for the doped Ceo and C70 [14] in 
the theoretical characterizations of the optical absorption experiments. 
Second, we will study with the model II [7] , 

H ^ -t J2 E(cW + h-C.) 

(i j):intralaycr 
(i,j):intcrlaycr 

+ Uj^^iA^iA^ (2) 

i 

where the sum of the second line is taken over pairs of the sites in neighboring layers with 
a cutoff for long distance as used in [15] of the calculation for multiwall carbon nanotubes; the 
function /5(r) is given by 

(3{r) = Aexp{-r/0 (3) 

where r is the distance between carbon atoms, A — 5.21t, and ( — 0.86A. The magnitude, 
(3(r — 3.40A) — OAt, is a typical value for the interlayer interaction strength in the tight 
binding model for A-B stacked graphite layers [16]: the exphcit value is about 0.35 - 0.39 eV, 
and t ~ 3eV gives the interaction strength about OAt. 

We investigate the systems which have open shell electronic structures when a nanographene 
layer is isolated. The effects of additional charges coming from functional side groups are 
simulated with introducing site potentials Eg [17,18] at edge sites in the model I or II. The 
form of the site potentials is 

H,i,^^ E,Y,Y.4,aCi,a, (4) 

where the sum of i is taken over the set of the positions of site potentials /. When Eg > 0, the 
electron number decreases from the average value at the site i, and the site potential means the 
electron attractive groups. When Eg < 0, the electron donative groups are simulated because 
of the increase of the electron number at the site potentials. Here, we take Eg = —2t, and one 
additional electron per layer is taken account. 



The finite size system, whose number of the stacking layers is similar to that of the samples 
with the nanometer size, is solved numerically with using the periodic boundary condition for 
the stacking direction, and we obtain two kinds of solutions. One of them is an antiferromagnetic 
solution, where the number of up spin electrons is larger than that of down spin electrons in 
the first layer, the number of down spin electrons is larger than that of the up spin electrons in 
the second layer, and so on. The other kind of solution is a nonmagnetic solution. The present 
author has discussed the antiferromagnetism in Ceo polymers, too [19]. The same technique 
(unrestricted Hartree-Fock approximation) used in ref. [6,7,19] is effective for the open shell 
electronic systems of this paper. The number of unit cells in the stacking direction is 10, and 
therefore there are 20 layers in the system used for the numerical treatment. 

In the model I, the parameters are changed within < ti < O.bt and < [/ < 4t. The 
realistic value of ti is estimated to be about O.lt at most [16], but we change this parameter 
for more extended regions in order to look at the behaviors of solutions in detail. In the 
antiferromagnetic solutions, the number of electrons is the same with the number of sites, and 
electronic states are half-filled. Because we assume the electron donative case Eg < 0, the 
number of electrons per layer is larger than the site number. Here, we take the electron number 
per layer = 25 and the site number per layer rigite = 24. All of the quantities of the energy 
dimension are reported using the unit t (~ 2.0 — 3.0 eV). 

In the model II, we investigate the continuous change between the A-B stacking [Fig. 1 
(a)] and the A-A stacking [Fig. 1 (b)]. We will move the first layer of Fig. 1 (b) to the upper 
direction. When the relative shift d — 0, the geometry is of Fig. 1 (b). As d becomes larger, 
the system changes from the A-A stacking to the A-B stacking. When d — a (a is the bond 
length in each layer), the system has the geometry of Fig. 1 (a). In increasing d, the edge sites 
feel weaker interactions from neighboring layers. Such the changes will give rise to variations 
in magnetic properties. The magnetism of the model II will be discussed with varying d and R 
{R is the interlayer distance) in section 4. 
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3. Results of model I 



Here, we consider the Hubbard-type model for systems which have open shell electronic struc- 
tures when a nanographene layer is isolated. The effects of additional charges coming from 
functional side groups are considered using the model I whose form is given by eq. (1). 

Figure 2 displays the absolute values of total magnetic moment per layer. Figures 2 (a), (b), 
and (c) are for the A-B stacking, and Fig. 2 (d) represents the case of the A-A stacking. The 
strengths of the Coulomb interaction U are shown in the figure captions. In Fig. 2 (a), the site 
potentials locate at the site D in the first layer [Fig. 1 (a)], and at the symmetrically equivalent 
site D' in the second layer. The site potentials exist at the sites E and E' in Fig. 2 (b), and they 
are present at the sites F and F' in Fig. 2 (c). The total magnetization is a decreasing function 
in these three figures. The decrease is faster in Figs. 2 (b) and (c) than in Fig. 2 (a). The sites 
E and F are neighboring to the site with the interaction ti, and thus the localized character of 
the magnetic moment can be affected easily in these cases. The decease of magnetization by 
the magnitude 30 — 40% with the water molecule attachment [3] may correspond to the case 
of Fig. 2 (b) or Fig. 2 (c). 

Next, we note that the magnetic ordering is not present for the A-A stacking case [shown in 
Fig. 1 (b)] whose numerical result is shown in Fig. 2 (d). This is owing to the fact that all the 
sites between the neighboring layers interact via the hopping integrals, and that the itinerant 
characters of the electrons are dominant strongly for the A-A stacking. The present results are 
in agreement with the calculations for the closed shell systems reported in the previous paper 
[6]. 

4. Results of model II 

In this section, the effects of site potentials and additional charges are investigated with the 
model II taking into account of the continuous relative shift between layers. The model eq. 
(2) has been used in [7], and is reviewed in the section 2 of this paper. 
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Figure 3 shows the absolute magnitude of the total magnetic moment per layer as a function 
of the relative shift d. The site potentials with the strength Eg = —2t are present at the sites D 
(E, and F) in the first layer, and at the sites D' (E', and F') in the second layer, in Figs. 3 (a) 
[(b), and (c)], respectively. Here, the interlayer distance R is fixed and the Coulomb strength 
U is varied in the series of the plots. The system is with the A-A stacking at c? = 0, and the 
stacking is of the A-B type at d — a. At d — and in the smaller region of d, there is not 
finite magnetization due to the strong interlayer hopping interactions near the A-A stackings. 
This qualitative property agrees with that of the model I and also with the calculation of the 
closed shell electron system [7]. However, there appear the finite magnetic moments with the 
antiferromagnetic sign alternation in the stacking direction for the region d > 0.5a and for U 
larger than a critical value. The onset of the magnetic moment at d — 0.5a might be due to 
the presence of a cutoff of the long distance interaction, which has been used in [15], also. The 
magnetic moment increases as U becomes larger as expected for the increase of the localized 
characters of electrons. In contrast, the magnetic moment is a weak decreasing function with 
respect to d. 

Figure 4 shows the same quantity as a function of d. In these plots, the Coulomb interactions 
are taken constant as a representative strength U — 1.2t, and the interlayer distance R is varied. 
The interlayer distance becomes smaller from the filled squares, through the open squares, to 
the filled circles. In the regions of d where the finite magnetization is present, the magnetic 
moment decreases strongly while the interlayer distance becomes shorter. The total hopping 
interactions become larger as R decreases. This effect enhances the itinerant characters of 
electrons in the stacking direction. Therefore, the magnetic moment becomes smaller with 
respect to the decrease of the interlayer distance. Such the qualitative behavior of the open 
shell electron system is in contrast with that of the closed shell electron system reported in 
the paper [7]. The present calculation is in qualitative agreement with that of the experiment 
which reports the decrease of the magnetization in the course of the adsorption of molecules 
[4,5]. In fact, the decrease of the interlayer distance about 0.4 A while the adsorption of water 
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molecules has been reported [20]. 

The magnetic moment decreases easily for the site potentials at E and E' [Figs. 3 (b) and 
4 (b)], and also for the site potentials at F and F's [Figs. 3 (c) and 4 (c)]. It does not vary so 
much in the case of the site potentials at D and D' [Figs. 3 (a) and 4 (a)]. Such the quantitative 
difference comes from the contrast whether the site potentials are far from the sites with strong 
interlayer hopping interactions or not. The difference of the effects of the site potentials at D 
(D') from those of the potentials at E (E') or F (F') agrees qualitatively with the result of the 
model I reported in the previous section. 

To conclude, the two calculations of models I and II of this paper agree with the experi- 
ments, qualitatively. The case of the site potential at E or F agrees with experiments, even 
quantitatively. We can explain the decrease of magnetism of 30 — 40 % in the process of ad- 
sorption of molecules [4,5]. Therefore, the open shell electronic structure due to the active side 
groups is a good candidate which could explain the exotic magnetisms. The exotic magnetism 
is related with the edge states which have large amplitudes at the zigzag edge atoms of the 
each graphene layer. The important roles of the edge states have been discussed intensively for 
the nanographite ribbon systems recently [21-23]. 

5. Summary 

Antiferromagnetism in stacked nanographite has been investigated with the Hubbard-type mod- 
els. The A-B stacking or the relative interlayer relation near the A-B stacking is favorable for 
the hexagonal nanographite with zigzag edges, in order that magnetism appears. We have 
also found that the open shell electronic structure can be an origin of the decreasing magnetic 
moment with adsorption of molecules. 
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Figure Captions 



Fig. 1. Stacked nanographite with zigzag edges. The bold and thin hnes show the first and 
second layers, respectively. The stacking is the A-B type in (a), and it is the simple A-A type in 
(b). There are 24 carbon atoms in one layer. The circles in (a) show sites where the interlayer 
distance between carbon atoms is shortest. The sites, A, B, and F, are edge sites of the first 
layer, and the sites. A', B', and F', are edge sites of the second layer. When there is a site 
potential at D (E or F) of the first layer, another site potential is present at D' (E' or F') of 
the second layer. 

Fig. 2. The absolute magnitude of the total magnetic moment per layer of the model I as a 
function of ti. The stacking is the AB-type in (a), (b), and (c), and it is the AA-type in (d). 
There is a site potential Es — —2t, (a) at the site D, (b) at the site E, and (c) at the site F. 
The site positions are displayed in Fig. 1 (a). In (a), the onsite interaction is varied within 0.6t 
(closed squares) < U < 1.8t (closed triangles). The interval of U between the series of the plots 
is AC/ = 0.3t. In (b), it is varied within l.Ot (closed squares) < U < 2.0t (closed triangles). 
The interval of U between the series of the plots is AU = 0.25t. In (c), it is varied within 0.6t 
(closed squares) <U < 1.8t (closed triangles). The interval of U between the series of the plots 
is AU = 0.3t. Finally, in (d), we display that there is no magnetization for the AA-stacking 
with a site potential at the D, E, or F. 

Fig. 3. The magnitude of the total magnetic moment per layer of the model II as a function of 
d and U. There is a site potential Eg — —It, (a) at the site D, (b) at the site E, and (c) at the 
site F. The site positions are displayed in Fig. 1 (a). The interlayer distance is R = 3.4A. The 
values of U are U — 0.4i (filled squares), 0.8i (open squares), 1.2t (filled circles), IM (open 
circles), 2.0t (filled triangles), respectively. 
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Fig. 4. The magnitude of the total magnetic moment per layer of the model II as a function of 
d and R. There is a site potential Eg = —2t, (a) at the site D, (b) at the site E, and (c) at the 
site F. The site positions are displayed in Fig. 1 (a). The strength of the Coulomb interaction 
is U — 1.2t. The values of R are R - 3.4A(filled squares), 3.2A(open squares), and 3.0A(filled 
circles), respectively. All the plots within < c? < 0.5a overlap, so only the filled squares are 
seen. 
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